Effects of Strain on Electronic Properties of Graphene 
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We present first-principles calculations of electronic properties of graphene under uniaxial and 
isotropic strains, respectively. The semi- metallic nature is shown to persist up to a very large 
uniaxial strain of 30% except a very narrow strain range where a tiny energy gap opens. As the 
uniaxial strain increases along a certain direction, the Fermi velocity parallel to it decreases quickly 
and vanishes eventually, whereas the Fermi velocity perpendicular to it increases by as much as 25%. 
Thus, the low energy properties with small uniaxial strains can be described by the generalized 
Weyl's equation while massless and massive electrons coexist with large ones. The work function is 
also predicted to increase substantially as both the uniaxial and isotropic strain increases. Hence, 
the homogeneous strain in graphene can be regarded as the effective electronic scalar potential. 
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Mechanical strain often gives rise to surprising effects 
on electronic properties of carbon nanomaterials 1 -"— . It 
can turn the metallic nanotube into semiconductor and 
vice versai~— . Along with the uniquely strong mechanical 
properties of the sp 2 - and sp 3 - bonded carbon materials 6 , 
the interplays between mechanical and electronic prop- 
erties may be useful in various applications 7 . A recent 
successful isolation of a new carbon allotrope^, graphene, 
offers a new opportunity to explore such interesting elec- 
tromechanical properties in two dimensions. 

At low energies, graphene at equilibrium has two lin- 
ear energy bands that intersect each other at the high 
symmetric points, K and K' , of the first Brillouin zone 
(BZ) and are isotropic with respect to the points^. With- 
out strains, the density of states vanishes linearly at the 
Fermi energy (Ep) or the Dirac point (Ed), exhibiting a 
semi-metallic nature. Thus, charge carriers are well de- 
scribed by the Dirac's equation for a (2+l)D free mass- 
less fermion2r— . Electron states here have another quan- 
tum number called a pseudospin which is either parallel 
or antiparallel to the wavevector of the electron and is 
of central importance to various novel phenomenal— . 
Mechanical strains can introduce new environments in 
studying such novel physics of graphene. 

Recently, several experiments have been performed to 
investigate the physical properties of graphene when its 
hexagonal lattice is stretched out of equilibrium^— . 
Strain can be induced on graphene either intentionally 
or naturally. The uniaxial strain can be induced by 
bending the substrates on which graphene is elongated 
without slippage 1 ^"—. Elastic responses are measured 
by pushing a tip of atomic force microscopes on sus- 
pended graphene 1 -. Graphene on top of SiO^ or SiC 
surface 2 ^ also experiences a moderate strain due to sur- 
face corrugations or lattice mismatch. Motivated by re- 
cent work s 14 ' 21 " — pointing to a remarkable stability of 
graphene with large strains, we have carried out first- 
principles calculations and theoretical analysis to explore 
the electronic structures of strained graphene and to un- 
derstand its low energy electronic properties. 



In this Rapid Communication, we show that no sizable 
energy gap opens in uniaxially strained graphene and the 
variation in energy bands strongly depends on the direc- 
tion of uniaxial strains. We also predict that the work 
function increases substantially as both the uniaxial and 
isotropic strain increases. When an uniaxial strain less 
than 26.2% is applied along the zigzag chain direction, 
the semi-metallicity is sustained. Beyond that, the sys- 
tem develops a small energy gap up to 45.5 meV at a 
strain of 26.5% and then close its gap quickly due to 
the downshift of the a* band to the Ep. This differs 
from conclusions of the previous literature s 14 ' 22 ' 23 . With 
uniaxial strain along the armchair chain direction, no en- 
ergy gap develops. Under uniaxial strain, the group ve- 
locities at the Ep are shown to be strong functions of 
the wavevectors so that the low energy properties with 
small uniaxial strains can be described by the generalized 
Weyl's equatio n 12 ' 13 ' 25 " — . With large uniaxial strains, 
quasiparticles become massive along the strain direction 
while ones in the perpendicular direction are still mass- 
less. 

Computations were carried out using the pseudopo- 
tential density functional method with a plane- wave basis 
set—. The exchange-correlation interactions were treated 
within the Perdew-Bcrkc-EnzelIio£2°» generalized gradient 
approximation. The cutoff energy for expansion of wave- 
functions and potentials was 400 eV and the Monkhorst- 
Pack k-point grid ofl2xl2xlis used for the atomic 
relaxation and of 60 x 60 x 1 for electronic structure calcu- 
lations. The atomic relaxation was carried out until the 
change in the total energy per one unit cell was smaller 
than 0.1 meV. The layer-to-layer distance between adja- 
cent graphene in the supercells is 15.0 A. 

Here, we consider graphene only under uniaxial and 
isotropic strains, respectively. For comparison, the elec- 
tronic structures of graphene under uniaxial strains along 
the two special directions are investigated. The effects of 
uniaxial strain along arbitrary directions and those of 
isotropic strains will also be discussed later. Following 
previous conventions^ 3 -, the uniaxial strain along the zig- 
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FIG. 1: (color online) (a) Hexagonal lattice of graphene. 
ai and 0,2 are the lattice vectors. With Z(A)-strain, 01 = 
(a x ,a y ) and <i2 = (— a x ,a y ). Si (i = 1,2,3) connects three 
nearest neighbors, (b) The first BZ with high symmetric 
points. Energy contours for graphene (c) without strain, (d) 
/1-strain of 20% and (e) Z-strain of 20%. The scale bar for 
contours is in unit of eV. The n and n* bands with various 
(f) ^4-strains and (g) ^-strains along the line of k y = in (b). 



zag chain direction [x-axis in Fig. 1(a)] in the honeycomb 
lattice is denoted by the Z-strain and one perpendicu- 
lar to this (y-axis) by the A-strain. From the fully re- 
laxed atomic geometries, the calculated Poisson's ratios 
for graphene as functions of the magnitude and direction 
of strains agree with the previous calculation s] 22 ' 24 . 

We find that if the magnitude of strain is less than 
26.2%, no gap opens with the Z-strain. Graphene with 
the j4-strain also has no energy gap up to a magnitude of 
30%. As shown in the energy contour from first-principles 
calculations, the Ed's coincide with the high symmetric 
K and K' (or R) points of the first BZ without strains 
(Fig. 1 (c)). With the dl-strain, the Ed's are off the sym- 
metric points and the two adjacent Ed's along the k y = 
line repel each other as the strain increases (Figs. 1(d) 
and (f)), agreeing with previous calculations 23 . Contrary 
to the cases with the dl-strain, the two adjacent Ed's 
with the Z-strain approach each other (Figs. 1(e) and 
(g)) and merge together eventually at strain of 26.2%. 

The mismatch of the Dirac points with the high sym- 
metric BZ points can be easily understood by one- 
orbital tight-binding approximation s 23 ' 25 ' 26 . In the elas- 
tic regime under the if-strain, the kinetic hopping inte- 
grals (t) between the nearest neighbors will depend on its 
connecting vectors, (5; (i = 1, 2, 3) such that t\ = £3 < t 2 
where U = t(Si) (i = 1,2,3) (Fig. 1(a)). Under the di- 
strain, t\ = £3 > t 2 . Considering the nearest- neighbor 
hoppings only, the Hamiltonian of graphene with Z(A)- 

strain can be written as W = — t 2 J2k £(k) c Ak CBk + c.c. 












27.0% 





FIG. 2: (color online) Calculated band structures of the 
strained graphene around the S point with a Z-strain of 
24.0%, 26.2%, 26.5% and 27.0% (from left to right panels), 
respectively. 



where £(k) = e k S2 (1 + cos(k x a x )), rj = t 1 /t 2 = 

ta/t2, k — (k x ,k y ) and ca(b)h is an annihilation oper- 
ator for an electron with momentum k on the sublat- 
tice A(B). The resulting energy dispersion is given by 
E k = ±t 2 \£(k)\. rj < 1 (77 > 1) for the Z(A)-strain. 
On the k y = line in the first BZ, the x-component 

of if -point is given by kx = [\ + ^ whereas the 

Dirac point with strains, i.e., the zero energy solution, 
£(k D ,0) = 0, is given by k D = ^cos _1 (-^). Hence, 
under the A(Z)-strain, kD kx as shown in Fig. 1. 

We find that the energy splitting between the a and 
a* bands at the S point is reduced when the Z-strain 
increases (Fig. 2) and one at the T point does with the 
A-strain (not shown here). The strain-induced small en- 
ergy gap is eventually closed due to downshift of the a* 
band at the Z-strain of 27% (Fig. 2). In very high strain 
regime, a single orbital tight-binding approximation fails 
to capture the downshift of er*-orbitals although it shows 
approximately similar variations of 7r-bands in the low 
and moderate strain regimes 2 ^,. It is noticeable that it 
(tt*) electrons along SR become massive but that those 
along ST are still massless after the gap closure (Fig. 
2). Anomalous area expansion, i.e., the negative Pois- 
son's ratio^L is found when the a* band touches the Ep 
at the Z-strain larger than 27% because the antibond- 
ing states are occupied (the unit cell area increases by 
35% under the Z-strain of 30%). However, at this point, 
graphene may not be stabl o 15 ' 22 . Hereafter, we will con- 
sider graphene with strains less than 26. 5%22. 

As uniaxial strain increases, the group velocity at the 
Ed increases or decreases substantially depending on the 
wavevectors (Fig. 3). We calculate the group velocities 
of electrons by differentiating the energy dispersion of 
conduction bands directly, i.e., i>k = j- [^^] E — Ef - The 
group velocity along the dl-strain (va3 in Fig. 3) de- 
creases as increasing strain while ones (vai and vaa) in 
direction perpendicular to strains increase. Up to the 
A-strain of 24%, va3 is reduced by almost 60% of the 
group velocity without strains (i>o) and vai and vaa in- 
crease linearly by 25%. We also find that vai differs 
vaa (opposite direction to the former) as shown in Fig. 
3(a). Along the specific direction 2 in insets of Fig. 3, 
vai — VZ2 — v o- Under the Z-strain, the similar behav- 
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FIG. 3: (color online) The group velocities (vAi) of n and ir* 
electrons with the A-strain (a) and vzi with the Z-strain along 
the direction i (= 1,2,3,4 in insets) in an unit of isotropic 
group velocity (vo) without strain. The angle between the 
direction A2 and A3 in inset of (a) is 52° and one between Z2 
and Z3 in (b) is 38°. 



iors occur (Fig. 3(b)). It is also noticeable that vz\ dif- 
fers vz4 (opposite direction to the former) and the both 
become zero at the strain of 26.2%. We note that the 
group velocity anisotropy under strains may lead to an 
anisotropy of resistance shown in a recent experiment 15 . 

The low energy properties of graphene with moderate 
strains as revealed by our first-principles calculations 
can be described well by the generalized Weyl's equa- 
tion^—. By expanding £(k) around (fco,0) up to the 
first order of small momentum q, £(q) = £,{kD+q x ,q y ) — 
{Arj 2 — l) x ^ 2 a x q x — ia y q y . The resulting Hamiltonian can 
be written as T-L ~ v x a x q x + v y a y q y where cr x r y ) are Pauli 
matrices, v x = t2a x (4rf — I) 1 / 2 and v y — t2a y . With 
the Z-strain, ti increases predominantly over a contrac- 
tion of ctjp^ so that v y increases. On the other hand, 
v x = t 2 a x (4rf-l) 1 / 2 = ha^-l/r] 2 ) 1 / 2 < V3ha x since 
77 < 1 with the Z-strain. Hence, v x decreases very quickly 
upon elongation of a x followed by reduction of t\ with the 
Z-strain. For the A-strain, the opposite situation occurs. 
We note that this Hamiltonian also describes the low en- 
ergy physics of graphene superlattic o 12 ' 13 and a-(BEDT- 
TTF)2l3^ r — respectively. Thus, like graphene superlat- 
tices^ 3 -, the pseudospin in uniaxially strained graphene is 
not in parallel or antiparallel to the wavevectors suggest- 
ing some interesting transport propertie a 12 ' 33 i 34 . 

Although the simple model described above can ex- 
plain the results from our first-principles calculations in 
general, we should point out that the next-nearest neigh- 
bor (nnn) hopping, t', plays an important role in the low 
energy properties^. As shown in Fig. 3, the group veloc- 
ities along +k x {vA(z)i) differ one along — k x {va{z)\) i m_ 
plying tilted anisotropic Dirac cones due to the nnn inter- 
actions 2 ^—. With the A(Z)-strain, t' also depends on its 
six connecting vectors, such that t' a = i'(±ai) = i'(±a 2 ) 
and t' = t'(±{a 1 - a 2 )). x = t' p /t' a < 1(> 1) for the 
Z(A)-strain. The effective Hamiltonian for the nnn in- 
teractions around {ko, 0) can be written as H' ~ v' x q x ao 
where v' x = a x t' a (1 — xM) (4— l/^ 2 ) 1 / 2 and ao is an iden- 
tity. The resulting energy dispersion can be expressed 
concisely as £ q = v(4>q)q where q = tan -1 (q y / 'q x ) , 



FIG. 4: (color online) Calculated density of states of graphene 
with (a) the A-strain and (b) Z-strain. 

9 = (Qx + ll) 1 ' 2 and «Oq) = K cos0 q ± (v 2 cos 2 q + 
^sin 2 ^) 1 / 2 ^^. So, v A1(4) = v x T < and v zl{i) = 
v x ± v' x as shown in Fig. 3. Hence, the Dirac cone is 
tilted in the k x direction regardless of the uniaxial strain 
direction. With the large strain (> 20%), v' x becomes 
negligible so that va(z)i — v A(z)i- And, the tilting effect 
disappears when the large uniaxial strain is applied. 

The density of states (D(E)) around the Ed increases 
gradually as the uniaxial strain increase while maintain- 
ing its linearity (Fig. 4). With large strains (> 20 %), 
D(E) shows an abrupt change depending on the direc- 
tion of strains. Figure 4 shows the calculated D(E) with 
\E— Ep\ < 0.6 eV from first-principles calculations. From 
the generalized Weyl's equation, D{E) — where 

7T V p, 

1/4 = i dct)Jv 2 {(j>^^. The strain-induced re- 
ductions in the averaged anisotropic group velocities (v 2 F ) 
will increase the slope of D(E) as shown in Fig. 4. The 
D(E) changes significantly when a* band is near the Ep 
with the Z-strain (Fig. 4(b)). With the large Z-strain, 
the merging of two Dirac points signals the van Hove 
singularities of the n and n* bands (24% case in Fig. 
4(b)). When the gap opens with the Z-strain of 26.5%, 
D(E < E F ) ~ y/E 25 and D(E > E F ) shows a steep 
enhancement due to the a* band. 

The work function in uniaxially strained graphene is 
predicted to increase substantially as the strain increases 
(Fig. 5). The calculated work function of graphene 
without strain is 4.5 eV agreeing with the previous the- 
oretical 3 ^ and experimental 3 - estimations. As the strain 
increases up to 12%, the work function increases lin- 
early by 0.3 eV regardless of the direction of strains as 
shown in Fig. 5. The work function rises up further 
to 5.2 eV as the A-strain reaches 26%. However, with 
larger Z-strains, the work function saturates to 4.8 eV. 
Hence the variations in the work function can also char- 
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FIG. 5: (color online) Calculated work functions ($) of 
graphene with the A- and Z-strain. 
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FIG. 6: (color online) (a) Calculated Fermi velocity variation 
under the /-strain in an unit of the Fermi velocity (vo) without 
strain, (b) Calculated work functions ($) the /-strain. We 
set the work function without strain to zero here. 

acterize the direction of the strain. Our calculated re- 
sults indicate that the controlled charge transfer between 
gaseous molecules and graphene can be realized by strain- 
ing graphene. We also anticipate that the strain affects 
the band lineup at the graphene-metal contact^. 

To study the effect of uniaxial strains in arbitrary 
directions, we study the band structure of graphene 
stretched along the direction rotated by 10.9° with re- 
spect to the x-axis in Fig. 1(a). We confirm that no en- 
ergy gap opens up to a strain of 30% (not shown here). 
The work function also increases as strain increases. Our 
ab initio calculations conclude that no energy gap opens 
under uniaxial strain less than 26% along any arbitrary 
direction. 

Finally, we calculate the variations of electronic prop- 
erties of graphene under the isotropic strain (/-strain). 



Because the /-strain maintains all crystal symmetries of 
graphene, the electronic structures show no significant 
changes unlike uniaxially strained cases. The Fermi ve- 
locity decreases linearly to 86% of «o as the /-strain in- 
creases up to 10% (Fig. 6 (a)). The work function of 
the system also increases linearly up to 0.64 eV as the 
/-strain reaches 10% (Fig. 6(b)). From the calculation 
results, it is shown that the uniform strain induces effec- 
tive vector— and electric scalar potential in graphene. 

In summary, from first-principles calculations, it is 
shown that strained graphene does not develop an energy 
gap and that the group velocities under uniaxial strain 
exhibit a strong anisotropy. We show that the gener- 
alized Weyl's equation is an appropriate model for uni- 
axially strained graphene that incorporates all assessed 
properties that go beyond the simple tight-binding ap- 
proximations. It is also shown that the work function 
of strained graphene increases substantially as strain in- 
creases. 

Note added.- After submission, we became aware of 
related work on similar systems from other group s 37 ' 38 . 
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